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Abstract. We have searched for near-infrared stellar 
counterparts of IRAS point sources in the Large Mag- 
ellanic Cloud (LMC), in J- and K-bands. This resulted 
in the detection of 21 counterparts, of which 19 are new 
discoveries. Using colour-magnitude and colour-colour di- 
agrams, we identify 13 Asymptotic Giant Branch (AGB) 
stars with thick circumstellar dust envelopes, 7 possible 
early post-AGB stars or stars recovering from a ther- 
mal pulse, and 1 red supergiant or foreground star. For 
10 of the IRAS targets we do not succeed in detect- 
ing and/or identifying a near- infrared counterpart. We 
serendipitously detect 14 other red sources, of which 2 
are known Long Period Variables, and a few galaxies. The 
near-infrared and optical colours of the galaxies may indi- 
cate considerable interstellar extinction through the LMC, 
as much as Ay ~ 2-4 mag. The relative number of AGB 
carbon stars over oxygen stars is shown to decrease as 
the luminosity increases. Yet amongst the faintest mass- 
losing AGB stars oxygen-rich stars still exist, which puts 
constraints on current convection theories that predict the 
occurrence of third dredge-up and Hot Bottom Burning. 
We investigate the nature of some LMC stars that have in- 
frared properties very similar to suspected Galactic post- 
AGB stars. 

Key words: Stars: mass loss - Stars: AGB and post- 
AGB - dust, extinction — Magellanic Clouds - Infrared: 
galaxies — Infrared: stars 



* based on observations collected at the European Southern 
Observatory, La Silla, Chile (proposal ESO 56.E-0681) 



1. Introduction 

Population studies of Asymptotic Giant Branch (AGB) 
stars in the Large Magellanic Cloud (LMC) have led to 
a number of interesting results, using the advantage of 
its known distance. The period-luminosity relation for 
Mira variables was deduced from photometric monitor- 
ing of Long Period Variables (LPVs) in the LMC (Feast 
et al. 1989). Objective prism survey studies of the LMC 
showed a clear lack of carbon stars amongst the most lu- 
minous AGB stars (—6 >Mboi > —7 mag; see Iben 1981 
for a review). These studies have been largely based on 
samples of optically visible AGB stars (e.g. Blanco et al. 
1980; Cohen et al. 1981; Westerlund et al. 1981). 

Stars near the tip of the AGB experience mass-loss 
rates as high as 10~ 5 M a yr _1 , giving rise to circumstel- 
lar envelopes (CSEs) in which dust forms, making them 
optically thick at optical wavelengths. Hence these AGB 
stars may only be observable at infrared (IR) wavelengths 
and were missed in earlier studies. Samples of candidate 
obscured AGB stars in the LMC have only recently be- 
come available (Reid 1991; Wood et al. 1992; Zijlstra et al. 
1996). The AGB nature of some of these stars still needs to 
be confirmed, and their (circum-)stellar properties are yet 
to be explored in detail. The small existing samples should 
be extended both in number and in (circum-) stellar pa- 
rameter space, to allow a population study of mass-losing 
AGB stars in the LMC. 

In paper I (Loup et al. 1997) we constructed an IRAS- 
selected sample of AGB and post-AGB star candidates. 
In paper II (Zijlstra et al. 1996) we confirmed the AGB 
nature for a subsample thereof, on the basis of their bolo- 
metric luminosities and their near-IR (NIR) and IRAS 
colours. We also showed that the chemical type — oxy- 
gen or carbon star — may be deduced from the position 
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Fig. 1. IRAS 12 magnitudes versus IRAS [12]-[25] colours 
for the sample of IRAS sources that we selected as targets 
for observation in the NIR (open circles), the targets that we 
actually observed (shaded circles), the AGB star candidates 
from paper II (dots), and the red supergiants from paper II 
(crosses) 



of a mass-losing AGB star in the K— [12] versus H-K di- 
agram. Considering the difficulty of obtaining spectra of 
these stars, the IR colour diagnostics are a powerful tool in 
determining the nature of the IR stars (e.g. van der Veen 
k Habing 1988; Guglielmo et al. 1993; Le Bertre 1993; Le 
Sidaner & Le Bertre 1994; Loup & Groenewegen 1994). 

Here, in paper III of our series on the study of obscured 
AGB stars in the Magellanic Clouds, we present the results 
of a NIR search for counterparts of IRAS point sources in 
the direction of the LMC, aiming at enlarging the sample 
of known obscured AGB stars in the LMC. After classify- 
ing the stars as being mass-losing AGB stars, or not, we 
briefly comment on their properties. We discuss not only 
the AGB stars, but also the nature of other positive iden- 
tifications, and serendipitous detections. Amongst them 
are stars that may be explained as being in the post- AGB 
phase, and we compare their properties with those of a 
sample of IRAS-selected post- AGB star candidates in the 
Milky Way. We also find a few galaxies, which may be used 
as probes to measure the interstellar extinction inside the 
LMC. 

2. Observations 

2.1. Sample 

For our search for NIR counterparts, we reduced the list 
of nearly 200 best-candidate sources from paper I to a list 



of 68 sources, excluding stars that were studied in pa- 
per II (except for one source). Our selection criteria are 
either IRAS flux densities S25 > £12, or 6*12 > 0.35 and 
S25 > 0.22 Jy. Limited by the available observing time, 
we finally observed 31 sources. They were selected mainly 
amongst the sources that were relatively bright at 12 /Ltm, 
had red [12]— [25] colours, and were not affected by cirrus 
at 60 fim. We adopt [12]= -2.51og(5i 2 /28.3) and [12]- 
[25]= -2.51og(S*i 2 /S , 25) x (6.73/ 28.3), with S 12 and 5 25 
the flux density in Jy in the IRAS 12 and 25 /im bands, 
respectively (IRAS Explanatory Supplement 1988). In 
the IRAS [12] versus [12]-[25] magnitude-colour diagram 
shown in Fig. 1, these 31 sources are represented by shaded 
circles, whereas the remaining 37 sources are plotted as 
open circles. Several sources coincide in this diagram, be- 
cause of the discrete values for the IRAS flux densities. 
The diagram shows that we mostly selected the brighter 
and slightly redder sources. We also plotted the AGB star 
candidates (dots) and red supergiants (crosses) from pa- 
per II. The present sample is fainter at 12 /xm and has 
a redder [12]-[25] colour on average than the AGB star 
candidates from paper II, that are themselves fainter at 
12 /im and have redder [12]-[25] colours than the super- 
giants. 

2.2. J and K-band imaging photometry 

We observed on the clear nights 1/2 and 2/3 January 1996, 
using the NIR camera IRAC2 at the ESO/MPI 2.2m tele- 
scope at La Silla, Chile. Lens C was chosen to cover a field 
of view of 133" x 133", with a pixel scale of 0.49". This 
combines the advantages of a large field of view and a rea- 
sonable sampling of the point spread function for doing 
photometry. For each IRAS point source field we did a se- 
quence of 12 images in the K-band filter, each consisting 
of ten 3-second exposures, shifted by 5" in right ascension 
with respect to the previous image. We repeated this pro- 
cedure backwards, using the J-band filter. The search is 
deepest in the ~ 75" x 130" centred at the IRAS point 
source. We compared the IRAC2 fields with the Digitized 
Sky Survey ( http:/ /archive.eso.org/dss/dss ), to determine 
the actual field centres. The absolute pointing of the tele- 
scope was found to be accurate to about 10", from a com- 
parison of the IRAC2 field centres and the IRAS positions 
used to point the telescope. 

We constructed frames that represent a good approxi- 
mation of the background emission, by median-averaging 
the shifted images within each sequence, rejecting the 
brightest pixels to avoid contamination by stars. The 
background-subtracted images were flat-field corrected us- 
ing flatfields obtained by taking an image of a screen in the 
dome, illuminated by a lamp dedicated for this purpose, 
and subtracting a similar image with the lamp off. The 
individual images were shifted to bring the position of the 
IRAS point source back in the centre, by integer number 
of pixels to ensure flux conservation. Finally these images 
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were added together. In this way, we could detect stars 
down to limiting magnitudes of J ~ 20 mag, and K ~ 18 
mag. The limiting magnitude varied from field to field by 
about a magnitude, mainly due to differences in the back- 
ground. 

The standard stars were observed by taking an im- 
age consisting of an average of thirty 0.6-second expo- 
sures, and another image shifted by 40". Their difference 
yields a background-subtracted frame, which we then flat- 
field. The standard stars were observed regularly during 
the night, covering the same air masses as the LMC. 

The J and K-band images were blinked to identify the 
J-K reddest sources. On these, circular aperture photome- 
try was done with an increasing aperture radius, to create 
a radial magnitude profile of the star. The same was done 
for the standard star. From the differential magnitude pro- 
file the magnitude of the program star can be estimated, as 
well as its accuracy. The standard stars used are HD38150 
(Jsaao = 8.210 mag, K S aao = 7.880 mag), SA94-702 (= 
GSC 00048-00918; J S aao = 9.246 mag, K SAA o = 8.289 
mag), and HD52467 (J S aao = 8.637 mag, K SA ao = 8.699 
mag). These magnitudes are in the SAAO system (cf. 
Carter 1990), and the IRAC2 magnitudes have been con- 
verted to the SAAO system using relations derived by Lid- 
man (1995) on a single test night: 

/Oirac2 \ ( 0.883 0.117 \ /Jsaao \ m 
^Ki RA c2 J ^-0.055 1.055 J \K SAAO ) 1 > 

The l-<7 errors are 0.012 and 0.015 on the coefficients for 
the JiraC2 and Kjrac2 magnitudes, respectively. The con- 
ditions during the observations were photometric. Extinc- 
tion corrections were found to be 0.05 mag air-mass -1 
in the J-band, and 0.03 mag air-mass -1 in the K-band 
on the first night, and twice as large on the second night. 
The relative humidity, dome temperature, and seeing were 
about 60—70%, 13°C, and 0.7—0.8" on the first night, 
and about 40—60%, 15— 12°C, and 0.7—1" on the sec- 
ond night. 

2.3. BVRi-band imaging photometry 

We observed on the clear night 24/25, and the partially 
cloudy night 30/31 December 1996, using the direct imag- 
ing camera at the Dutch 0.9m telescope at La Silla, Chile. 

On the first night, we imaged two 3.8' x 3.8' fields cen- 
tred on the galaxies that we detected in the NIR close to 
the IRAS point sources LL-LMC0603 and LI-LMC1818 
(see below). The night was photometric, although the 
Moon was close to opposition. We took six images in the 
Bessel B-band, and three images in the Bessel R-band, 
all 200 seconds integration time per frame. Stellar images 
on these frames had FWHM of typically 1.1"-1.4". The 
images were calibrated by stars in the Landolt (1992) stan- 
dard star field SA98 that were observed close in time and 
airmass to the galaxy fields. Hence the B and R mag- 
nitudes are on the photometric systems of Johnson and 



Kron-Cousins, respectively. We obtained integrated mag- 
nitudes for the galaxies and the few rcdmost stars in the 
fields, using aperture photometry and applying aperture 
corrections derived from the much less crowded standard 
star fields. 

On the second night, we measured the bright stellar 
counterpart (see below) of LI-LMC1821 in the Bessel B, 
V, and R-bands, and the Gunn i-band. The night was 
cloudy, and useful photometry could only be obtained by 
rapidly switching back and forth between the star and 771 
Dor (= HR2194), an AO Main Sequence star of V= 5.7 and 
B-V= —0.03 mag at 1.5 degree separation, until stable 
and well-defined magnitudes could be determined. The B, 
V, and R magnitudes that we obtained for LI-LMC1821 
are on the Johnson photometric system, whereas the i 
magnitude is on the Gunn photometric system. We esti- 
mated the Gunn i-band magnitude for r/i Dor to be i= 5.74 
mag. Integration times were 0.3, 0.2, 0.2, and 0.5 seconds 
for BVRi respectively. To avoid problems with shutter de- 
lay times for these short integration times, we used the 
same integration times for both stars, and they were al- 
ways positioned at the same place on the CCD. 

3. Results 

3.1. Identification with IRAS sources 

The positive and tentative identifications, non-detections, 
serendipitous detections, and galaxies are listed in Table 
1, together with the IRAS flux densities and NIR mag- 
nitudes, where applicable. The numbers of the stars refer 
to the LI-LMC numbers in the Schwering & Israel (1989, 
1990) catalogue of IRAS point sources in the direction of 
the LMC. A suffix b, c, or d was added in the case of field 
stars or galaxies, neither of which were considered to be 
the counterpart of the quoted LI-LMC source. Positive 
and tentative identifications were found by first selecting 
the object in the field of view with the reddest J-K colour. 
Stars with J-K> 2 mag were considered to be the counter- 
part of the IRAS source, as no two such red objects ever 
appeared in the same field of view. In the absence of ob- 
jects with J-K> 2 mag, objects were considered to be the 
counterpart of the IRAS source, if they either had J-K> 1 
mag (i.e. cool and/or reddened stars, see e.g. Bell 1992) or 
were very bright (e.g. LI-LMC1821), and not too distant 
from the position of the IRAS source (i.e. within ~ 20"). 
The counterparts with IR colours similar to those of mass- 
losing AGB stars have been enumerated under positive 
identifications, as this can be considered a confirmation 
of their identification with the IRAS source. The other 
counterparts have been enumerated under tentative iden- 
tifications, and probably are stars with detached CSEs. 

LI-LMC0109 was first detected in November 1995, us- 
ing the NIR photometer at SAAO. IRAS05003-6712 was 
already known to be a mass-losing AGB star in the LMC 
(paper II). LI-LMC0530 is identified with SHV05 10004- 
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Table 1. Names (see text), NIR positions, separations A of IRAS and NIR source (in arcsec), J- and K-band magni- 
tudes and J-K colours, IRAS 12- and 25 /im fluxes (in Jy) and [12]-[25] colours for the objects in our sample. We adopt 
[12]-[25]= -2.51og(S , 12 /5 25 ) x (6.73/28.3), with S 12 and S 25 the flux density in Jy in the IRAS 12 and 25 ^m bands, respec- 
tively (IRAS Explanatory Supplement 1988). 1-a error estimates are given for the NIR photometry. 

LI-LMC IRAS RA(2000) Dec(2000) A J aj K q K J-K a 3 - K S12 S25 [12]- [25] 
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tentative identifications: post- AGB or thermal pulse stars in LMC 
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Fig. 2. Cumulative distribution over separation of the IRAC2 
and IRAS positions, of positive and tentative IRAC2 identi- 
fications with H-K> 3 (bold solid) and H-K< 3 (thin solid), 
serendipitously detected NIR point sources (dotted) and galax- 
ies (dashed) 



692755, an LPV with a period of 169 days (Hughes & 
Wood 1990, who measured J= 12.63 mag and K= 11.51 
mag). We identify LI-LMC1721b with SHV0547489- 
704450, an LPV with a period of 264 days (Hughes & 
Wood 1990, who measured J= 12.86 mag and K= 12.69 
mag): LI-LMC1721 is not associated with this LPV, be- 
cause the NIR counterpart that we detect for this IRAS 
source is much redder than the LPV. The same is true for 
LI-LMC1721c, whereas LI-LMC1721d is too faint to be 
the LPV. We identify LI-LMC0937b with SHV0523536- 
700128, an LPV with a period of 381 days (Hughes & 
Wood 1990, who measured J= 10.84 mag and K= 9.48 
mag). The IRAS source LI-LMC1759 lies at the edge of a 
small open cluster of perhaps a dozen stars within 0.5' di- 
ameter, and may therefore not be related to any individual 
star. The galaxies are discussed below. 



3.2. Serendipitous detections 

We have calculated the separations of the IRAS point 
sources and the NIR sources detected with IRAC2 (Fig. 
2: accumulated starting at infinity, and normalised to the 
total number of sources) . The IRAS position is generally 
accurate to 5-15 arcsec, depending mainly on the 12 /um 
flux level, but can be off by more than half an arc-minute 
(paper II). 

The reddest sources with J-K> 3 (bold line) arc un- 
doubtedly the IRAS counterparts. The combined area of 



our 31 fields covers ~ 10~ 3 square degrees. The area of the 
LMC covered by the selection in paper I measures ~ 10 2 
square degrees. Only if there were of the order of 10 5 such 
NIR-red sources would one expect to have detected one 
serendipitously. There are of the order of 10 2 IRAS de- 
tected (post-) AGB star candidates (paper I). Hence the 
population of stars with J-K> 3 that are not detected by 
IRAS would have to be at least 10 3 as many as the IRAS- 
detected stars in the LMC. The number of planetary neb- 
ulae (PNc) in the LMC is estimated at 1100 (Pottasch 
1984). The PN lifetimes of ~ 10 4 years are very similar 
to the lifetimes of the mass-losing AGB phase (less than 
10 5 years). Assuming that all mass- losing AGB stars will 
eventually form PNe, we estimate the number of mass- 
losing AGB stars in the LMC to be 10 4 at most, arguing 
against the possibility of a population of 10 5 LMC stars 
with J~K> 3. 

The bluer sources (J-K< 3) that we consider as posi- 
tive or tentative identifications (thin line in Fig. 2) could 
be contaminated by NIR stars that are not associated with 
the IRAS source. But these cannot be many, as seen from 
the similarity of the separation distribution to that of the 
reddest sources. The serendipitously detected NIR sources 
(field stars) and the galaxies have very similar separation 
distributions, much broader than the separation distribu- 
tions of the positive identifications. This is indicative of 
their not being the IRAS counterpart. 

A field of view of 130" x 75" (the deep area only) im- 
plies an expected mean separation of 40.6" for a serendipi- 
tous detection. The mean separation for our serendipitous 
detections is 37", and 38" for the galaxies. We suspect 
that the detection probability near the edges is somewhat 
lower. Essentially all our detections are situated in the 
deep area. The mean separation is 18" for the reddest, 
and 23" for the bluer positive and tentative identifications. 
Replacing 2 or 3 of the reddest stars with separations of 
18" by the same number of stars but with separations of 
37.5", we generate a separation distribution with a mean 
of <~ 23". Hence we expect there may be 2 or 3 serendip- 
itous detections in the group of positive and, more likely, 
tentative identifications with J-K< 3. This is too small to 
affect the conclusions that we reach in the present study. 

3.3. Galaxies as a probe of the interstellar extinction in- 
side the LMC 

We discovered a few galaxies, which are probably not re- 
lated to the IRAS point sources. Since the galaxies are 
located behind the LMC, they could, in principle, be used 
as probes of the interstellar reddening through the LMC. 
If the J-K colour excess of ~ 0.7 mag is due to extinc- 
tion by dust inside of the LMC, a visual extinction of 
Ay ~ 4 mag is indicated; the extinction through the en- 
tire LMC would vary between Ay ~ and Ay > 4 mag. 
This could have severe consequences for the observation of 
stars inside of the LMC. Although Oestreicher et al. (1995) 
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Fig. 3. J-band (left) and K-band (right) images of the edge-on spiral galaxy LI-LMC1818d. Coordinates are in arc-seconds, 
relative to the nucleus of the galaxy. Contour levels are between 2.2 x 1CT 28 — 1.1 x 1CP 26 erg cm -2 s _1 Hz -1 per square 
arc-second, with intervals of 2.2 x 1CP 28 erg cm~ 2 s _1 Hz -1 per square arc-second, the same for both images. Note that East 
is down and North is to the left 
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Fig. 4. J-band (left) and K-band (right) images of the face-on spiral galaxy LI-LMC0603c (A), interacting with a second (B), 
and possibly a third galaxy (C). Coordinates are in arc-seconds, relative to the nucleus of galaxy A. Contour levels are between 
2.6 x 10" 28 — 6.4 x 10" 27 erg cm s Hz per square arc-second, with intervals of 1.3 x 10 erg cm s Hz per square 
arc-second, the same for both images. Note that East is down and North is to the left 



showed the foreground reddening towards the LMC to be 
only Ay ~ 0.18 on average, stars inside the LMC could 
suffer from visual extinction of a few magnitudes. 

Are the measured colours of the galaxies intrinsic to 
them, or have they been severely affected by interstellar 
reddening through the LMC? The intrinsic colours of a 
galaxy depend on its kind. The magnitudes for the galaxies 
approximately represent the total integrated light. There 
is no sign of steep colour gradients, although the bulges 
or nuclei appear somewhat redder. LI-LMC1818d (Fig. 

3) is a bright edge-on spiral galaxy. LI-LMC0603c (Fig. 

4) is a face-on spiral galaxy (A) interacting with another 
galaxy (B), and possibly with a third galaxy (C). All other 
objects in the images of LI-LMC1818d and LI-LMC0603c 
are unresolved. LI-LMC1759d is probably a spiral galaxy 
seen under a small inclination angle, but its position close 
to the edge of the field severely degraded the quality of its 
image. LI-LMC1803b is small and barely resolved. 

Normal galaxies are confined to a colour J-K ~ 
0.8—0.9 mag (e.g. Glass 1984; Silva 1996). Only LI- 
LMC1759d has a colour consistent with a normal galaxy. 



LI-LMC1818d may suffer from extinction by dust in its 
disk because of its edge-on orientation, but the red colour 
of the near-face-on LI-LMC0603c cannot be explained eas- 
ily by extinction inside of a normal galaxy. 

Galaxies with anomalous IR colours include emission 
line galaxies. Whitelock (1985b) presented NIR data on a 
sample of thirteen IRAS galaxies, probably H n galaxies. 
Their mean J-K colour was 1.25 mag, with a standard 
deviation of only 0.13 mag. But the J-K colours of our 
reddened galaxies are still consistent with the reddest H II 
galaxies, while both Seyfert 1 and 2 galaxies can reach J- 
K colours in excess of 2 mag (Glass & Moorwood 1985; 
Almudena et al. 1996; Kotilainen & Ward 1994). Active 
galaxies usually have strong colour differences between the 
nucleus and the rest of the galaxy, something which we do 
not see in our galaxies. 

There is less doubt about the question whether our 
galaxies could be the counterparts of the IRAS sources 
in the field. The large distances between the galaxies and 
the IRAS positions and/or the presence of a NIR-redder 
point source in the field already suggest that the detec- 
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tion of the galaxies was a mere coincidence. Also, the 
mid-IR colours of galaxies detected by IRAS are different 
from those of our sources. The IRAS galaxies discussed 
by Whitelock (1985b) have IRAS colours typical of cold 
dust, with S60/S25 ratios ~ 8 ± 1. This would have lead 
to 60 /xm flux densities for our IRAS sources of ~ 3 — 6 
Jy whereas the measurements indicate they can only be 
~ 1 Jy at most. Whitelock's galaxies have K magnitudes 
of ~ 11 — 12 mag, which is considerably brighter than our 
galaxies, whereas their IRAS 25 /um flux densities are com- 
parable. Hence we cannot exclude the possibility that our 
reddened galaxies have mid-IR excess emission below the 
sensitivity of IRAS, but they are not the counterparts of 
the IRAS sources. 

The B and R-band photometry for the red sources in 
the fields of LI-LMC0603 and LI-LMC1818 is presented 
in Table 2. The stars that we identified as the NIR coun- 
terparts of the IRAS point sources LI-LMC0603 and LI- 
LMC1818 were not detected down to - 22 nd mag in the 
R-band. The other NIR-rcd objects do indeed have op- 
tically red counterparts, both the stars and the galaxies. 
Moreover, we discovered two very red stars in the Dutch 
telescope field of LI-LMC0603c, and a red galaxy in the 
Dutch telescope field of LI-LMC1818d, all three which 
are outside the corresponding IRAC2 fields. The B-R 
colour of ~ 5 mag for the two red stars LI-LMC0603d 
and LI-LMC0603e implies severe inter- or circumstellar 
extinction (Whitelock et al. 1996). The red B-R colour 
of galaxy LI-LMC1818e is remarkably similar to that of 
the galaxy LI-LMC1818d. Active galaxies are often opti- 
cally blue (Veron-Cetty & Veron 1996 adopt B-R = 0.57). 
They can be optically red, though, as is the case for IRAS 
galaxies (P. Veron, private communication; see also Due 
et al. 1997). 

We have taken the Lauberts-Valentijn ESO catalogue 
of galaxies (Lauberts & Valentijn 1989) to investigate the 
variation in galaxy B-R colours across the sky in the vicin- 
ity of the LMC. We had to sample in four square degree 
bins in order to obtain useful statistics. We identified a 
region of 64 square degrees centred at 4 h Right Ascension 
and —62° Declination, that was relatively well populated 
(148 galaxies), and that appeared to be representative for 
the colours of galaxies unaffected by extinction through 
the Magellanic Clouds system. This region yielded a mean 
B-R= 1.09 ± 0.17 mag, with minimum and maximum B- 
R of 0.70 and 1.38 mag respectively. We can compare this 
to the bins containing our galaxies LI-LMC0603c and LI- 
LMC1818d. Each of these bins contained one galaxy from 
the above mentioned catalogue. Although this is statisti- 
cally very poor, it is interesting to note that the galaxy 
near LI-LMC0603c has B-R = 1.74 mag, and the galaxy 
near LI-LMC1818d has B-R =1.44 mag, both consider- 
ably redder than the estimated typical colour. Comparing 
the B-R colours of the galaxies LI-LMC0603c and LI- 
LMC1818d with the canonical value of B-R- 1.09, we 
arrive at extinctions of Ay — 2.2 mag. 
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If we compare the positions of our galaxies with the 
dust column density maps by Schwcring (1989), we notice 
that the dust distribution in the LMC is rather patchy 
as compared to the SMC. Two of our reddened galaxies 
are situated at the East side. Although they lie outside of 
the area covered by the dust map, it is at the East side 
that there appears to be a massive dust complex. Some 
of the serendipitously discovered red stars may thus be 
normal stars, situated behind a local, but large amount of 
dust in the LMC. Alternatively, they may be the nuclei of 
normal galaxies as seen through the LMC. In that case, 
the dust causing the reddening of the galaxy colours may 
be situated behind the main body of the LMC, leaving the 
colours of the stars in the LMC unaffected. 

The colours and morphologies of the galaxies, the fact 
that we detect only a few galaxies of which most are red, 
and the presence of other red objects in their projected 
vicinities all provide circumstantial evidence for severe in- 
terstellar extinction inside or behind the LMC. However, 
it is based on very little data, and should be confirmed 
by the systematic study of the colours of galaxies seen 
through the LMC. 

4. The nature of the IRAS counterparts 

In this section, we classify the IRAS counterparts as mass- 
losing AGB stars that are oxygen rich (O), carbon rich 
(C), or that are not distinguishable between oxygen stars 
and carbon stars (OC). We argue that some of the stars 
may be post-AGB or thermal pulse stars (post/TP). The 
arguments used in classifying the stars are explained be- 
low, and the results are tabulated in Table 4. 

4-1. Chemical classification from colours and magnitudes 
4.1.1. Galactic comparison sample 

In paper II we show how in a K-[12] versus H-K diagram 
carbon stars are distinguished from oxygen stars. A simi- 
lar diagram can be made using J-K. In Fig. 5 we present 
the K-[12] versus J-K diagnostic diagram, for the carbon 
stars (Fig. 5a) and the oxygen stars (Fig. 5b) from the 
Galactic sample described in paper II. It is clear that in 
the Milky Way the K-[12] versus J-K diagram can also 
be used to separate carbon stars from oxygen stars. The 
Galactic carbon star sequence (dotted straight line) is de- 
scribed by the empirical relation 

J-K=|(1 + (K-[12])) (2) 

The Galactic oxygen star sequence (solid curved line) ap- 
proximately satisfies an infinite series of the form 

00 

J-K= 13 x ^10~ 3 "(K- [12]) 3 ™ (3) 

n=0 

There is a hint a secondary sequence of oxygen stars, in- 
dicated by the dashed line in Fig. 5b. Later we discuss the 



8 Jacco Th. van Loon et al.: Obscured AGB stars in the Magellanic Clouds III 

Table 2. Positions, R-band magnitudes, and B-R colours for the red sources in the Dutch telescope fields of LI-LMC0603 and 
LI-LMC1818. 



LI-LMC 


RA(2000) 


Dec(2000) 


R 


B-R 


remarks 


0603b 


05 12 04.5 


-70 33 14 


18.86 ± 0.04 


3.4 ± 0.9 


NIR-red star 


0603c 


05 12 05.4 


7n oo CiA 

- 10 32 04 


i t _i_ n no 
1 ( .50 ± 0.02 


o on I n i n 

2.30 ± U.12 


face-on NIR galaxy 


0603d 


05 11 51.5 


-70 33 32 


15.45 ± 0.02 


5.06 ± 0.17 


red star outside IRAC2 field 


0603e 


05 12 23.4 


-70 33 05 


15.26 ± 0.02 


5.01 ± 0.11 


red star outside IRAC2 field 


1818b 


06 02 11.8 


-72 26 41 


19.49 ± 0.06 


3.7 ± 1.5 


NIR-red star 


1818c 


06 02 14.5 


-72 27 43 


19.62 ± 0.08 


2.7 ± 0.4 


NIR-red star 


1818d 


06 02 12.8 


-72 27 36 


17.77 ± 0.03 


2.44 ± 0.24 


edge-on NIR galaxy 


1818e 


06 01 54.0 


-72 27 04 


16.78 ± 0.02 


2.47 ± 0.10 


galaxy outside IRAC2 field 



1 



cv 



(a) 



m 1 — — i r 




carbon stars 



1 



(b) 



n 1 1 — — i r 



4 

J-K 




Fig. 5. K-[12] versus J-K diagram for the carbon stars (a) and oxygen stars (b) of the Galactic sample of Guglielmo et al. (1993). 
We adopt [12]= — 2.51og (S12/28.3), with S12 the flux density in Jy in the IRAS 12 /xm band (IRAS Explanatory Supplement 
1988). From these data, we derived the carbon star (dotted) and oxygen star (solid) sequences. The dashed line indicates the 
possible existence of a secondary sequence for the oxygen stars 



stars at this side of the principal sequence in the K-[12] 
versus J-K diagram in more detail. From now, we assume 
that there be no substantial differences between the K- 
[12] versus J-K diagrams for mass-losing AGB stars in 
the Milky Way and in the LMC. 

4.1.2. LMC sample 

The K-[12] versus J-K diagram for the stars of the present 
sample in the direction of the LMC is given in Fig. 6, to- 
gether with a diagram of the K-band magnitudes versus 
J-K. The sequences we have drawn in the latter diagram 
are estimated from the positions of AGB stars and red su- 
pergiants in the LMC as compiled by Loup & Groenewe- 
gen (1994). LMC red supergiants and foreground stars are 



predominantly occupying a linear sequence, approximated 
in Fig. 6a by: 

K = —2.5 + 10( J — K) (4) 

Mass-losing AGB stars in the LMC follow the curved se- 
quence towards red J-K colours, approximated in Fig. 6a: 

K = 10.3 + 0.4 ((J- K) -if (5) 

Using the K-[12] versus J-K diagram, we can classify 
several stars as being mass-losing AGB stars either on the 
oxygen star sequence (solid), or on the carbon star se- 
quence (dotted) , with the carbon stars to be found exclu- 
sively amongst the optically thickest sources. But there are 
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J-K 

Fig. 6. K-band magnitudes (a) and K-[12] colours (b) ver- 
sus J-K colours for the identified NIR counterparts of the 
IRAS sources (squares), the serendipitously detected stars (tri- 
angles), and the detected galaxies (circles), with la error bars. 
We adopt [12]= -2.5 log (Si 2 /28.3), with S 12 the flux density 
in Jy in the IRAS 12 fim band (IRAS Explanatory Supplement 
1988). The dotted lines in the upper panel represent the AGB 
stars in the LMC from Loup & Groenewegen (1994). In the 
lower panel the Galactic AGB carbon star (dotted) and oxy- 
gen star (solid and dashed) sequences from Fig. 5 are indicated. 
We also plotted arrows corresponding to an extinction in the 
visual of Av — 10 



several stars that have too large a 12 /jm excess for their 
J-K colour, or alternatively are too blue in J K for their 
12 fim excess to be on the AGB. Of these, LI-LMC1821 
is far too bright in the K-band to be an AGB star in the 
LMC, and it is probably a foreground star. Another pe- 
culiar source is LI-LMC0530, which lies on the Loup & 
Groenewegen (1994) sequence at a blue J-K colour. We 
will come back to this source later. The remaining six out- 
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Fig. 7. J-K (a) and K-[12] (b) colours, and IRAS 
12 (jm magnitudes (c) versus IRAS [12]-[25] colours for 
the identified NIR counterparts of the IRAS sources 
(squares), and the non-identifications (triangles), with 
1-cr error bars. We adopt [12]= -2.5 log (5i 2 /28.3) and 
[12]-[25]= -2.51og(S , i 2 /S , 2 5 ) x (6.73/28.3), with S12 and S 25 
the flux density in Jy in the IRAS 12 and 25 /im bands, re- 
spectively (IRAS Explanatory Supplement 1988). In the upper 
panel, the Galactic AGB carbon star (dottedO and oxygen star 
(solid) sequences from Le Bertre (1993) and Le Sidaner & Le 
Bertre (1994) are indicated 



liers are all faint in the K-band, and five of them are redder 
than J-K= 2. Up to three of the stars with lower limits to 
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their J-K colours might be similar to the group of six out- 
liers as well. The field stars are distributed over a larger 
range of K-band magnitudes, and like the galaxies they 
are all bluer than J-K= 2. This suggests that the outliers 
and the field stars are of different nature. 

Another way of separating AGB carbon stars from 
AGB oxygen stars is the J-K versus [12]-[25] colour- 
colour diagram. For this purpose, we have estimated the 
carbon star and oxygen star sequences from Le Bertre 
(1993) and Le Sidaner & Le Bertre (1994). The average 
carbon star sequence may be approximated by: 

[12] - [25] = -0.1 + 0.16(J-K) (6) 

The average oxygen star sequence may be approximated 
by: 

[12] - [25] = 0.55 + 0.21(J - K) (7) 

The J-K versus [12]-[25] diagram for the stars of our 
sample is given in Fig. 7, together with diagrams of the 
K-[12] colours and the 12 /xm magnitudes versus [12]-[25]. 
For the errors on the IRAS flux densities we have adopted 
formal values of la = 0.03 Jy. The J-K versus [12] -[25] di- 
agram does not work well for our sample mainly because 
the lower limits to the J-K colours of several stars allow 
these stars to lie either on the oxygen star sequence (solid) , 
or on the carbon star sequence (dotted). The outliers of 
the K-[12] versus J-K diagram have relatively blue J-K 
colours. Alternatively they may be characterised by rela- 
tively cool dust envelopes because of their relatively red 
[12]-[25] colours. From the [12] versus [12]-[25] diagram 
it appears that the IRAS sources for which no counter- 
part was found (triangles) may have cooler dust envelopes 
and/or may experience less severe mass loss, but these in- 
dications are very marginal (see also paper II). We also 
note that if the counterpart would not have been partic- 
ularly red or bright as compared to the other stars in the 
field, we may have failed to recognise the counterpart as 
such. This may be the case for PNe, which can have blue 
J-K colours between and 0.5 mag due to strong He I line 
emission at 1.083 ^m (Whitelock 1985a). 

4-2. Post- AGB star candidates 

It is interesting to compare the identified NIR counter- 
parts of the IRAS sources as found in the present study, 
with the sample from van der Veen et al. (1989: VHG-89). 
The latter sample consists of a compilation of Galactic 
objects thought to be in the transition from the AGB to 
the planetary nebula phase. Their distances and absolute 
magnitudes are not known, so that we can use this sam- 
ple in the colour-colour diagrams but not in the colour 
magnitude diagrams. 

The J-K, K-[12], and [12]-[25] colours arc plotted ver- 
sus each other in Fig. 8 (circles for the present sample, and 
dots for the VHG-89 sample). The VHG-89 stars have 



cool dust envelopes, and have larger 12 ^m excesses for 
their J-K colours than would be reconcilable with mass- 
losing AGB stars. The six outliers of the present sample 
roughly overlap with the VHG-89 sample. We classify the 
six outliers therefore tentatively as post- AGB stars. They 
would be related to the VHG-89 post- AGB stars with rel- 
atively warm dust envelopes, and either have relatively 
large 12 /im excesses or be optically thick in the K-band. 
We cannot exclude that one or more of the LMC post- 
AGB star candidates are actually AGB stars that have 
recently experienced a thermal pulse (see section Discus- 
sion). 

4.3. The stellar counterpart of LI-LMC1821 

The IRAS point source LI-LMC1821 is identified with a 
bright NIR star. Our BVRi photometry is presented to- 
gether with the NIR and IRAS photometry in Table 3. 
The B-V and V-i colours indicate a non-reddened, early- 
M type spectrum (Iyengar & Parthasarathy 1997). We 
plot the spectral energy distribution in Fig. 9 (squares), 
together with an arbitrarily scaled 3000 K blackbody. We 
also plotted mean fluxes of the M0 (open circles) and M5 
(solid circles) giant from Fluks et al. (1994), scaled to 
a distance of 1.5 kpc. LI-LMC1821 may be a late-K or 
carly-M giant at a similar distance, and a member of the 
Galactic halo or thick disk population (Robin et al. 1996). 
Alternatively, it may be a very luminous red supergiant 
in the LMC, with a bolometric luminosity of Mboi ~ —9.7 
mag. The small IR excess indicates a thin CSE as a result 
of modest mass loss. The IRAS colours are consistent with 
a detached CSE (van der Veen & Habing 1988), so that 
it may be a Galactic post-AGB star, but as the source is 
not detected at 60 /im this is not conclusive. It is surpris- 
ing that it has not been listed in any stellar catalogue or 
survey, despite that it is bright at optical wavelengths. It 
would be interesting to obtain a radial velocity measure- 
ment discriminating between membership of the LMC or 
the Milky Way. 

5. Bolometric luminosities 

We calculate bolometric luminosities for the stars that 
we identified as the NIR counterparts of the IRAS point 
sources. The method for estimating the unobserved part 
of the spectral energy distribution is based on the average 
NIR colours for obscured AGB stars, and will be described 
in a forthcoming paper. The method for obtaining the in- 
tegrated luminosity under the spectral energy distribution 
is described in paper II and originally in Whitelock et al. 
(1994). We adopted a distance modulus to the LMC of 
(m-M)o = 18.47 (Feast & Walker 1987). The results are 
tabulated in Table 4. 

The cumulative luminosity distributions of the differ- 
ent classes of stars are presented in Fig. 10a. In principle 
we could derive the luminosity distribution by differenti- 
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Fig. 8. K-[12] colours versus [12]-[25] (a) and J-K (b) colours, and J-K colours versus [12]-[25] colours (c), for the 
identified NIR counterparts of the IRAS sources in the present sample (circles, and squares for the possible post-AGB 
or Thermal Pulse stars) and the sample of post-AGB stars of VHG-89 (dots). We adopt [12]= — 2.5 log (£12/28. 3) and 
[12]-[25]= -2.51og(S , i 2 /5 , 2 5 ) x (6.73/28.3), with S12 and S25 the flux density in Jy in the IRAS 12 and 25 ^m bands, re- 
spectively (IRAS Explanatory Supplement 1988). The Galactic AGB carbon star (dotted) and oxygen star (solid and dashed) 
sequences are as in Fig. 6 and Fig. 7 



ation of the cumulative distribution function. In practice 
the errors would be huge because of the small number of 
stars. Instead, we apply Gaussian broadening by replacing 
each star with absolute bolometric magnitude Mbol,* by 

n = n norm x exp -(Mbol - M bo i,*) 2 (8) 

where n norm = 0.564 to normalise to unity per star. The 
FWHM of the star has thereby become 0.83 magnitude. 
Then we take the sum of the luminosity-broadened stars. 
A scientific motivation for applying the broadening to 
mass-losing AGB stars is that they are LPVs with bolo- 



metric amplitudes of about one magnitude but only ob- 
served at a single NIR epoch. The luminosity distributions 
are presented in Fig. 10b (line types are the same as in 
Fig. 10a). They represent the distribution over luminosity 
of the IRAS-detected AGB stars in the LMC, and they 
are likely to be incomplete, especially at the faint end. We 
assume that this does not seriously affect the relative dis- 
tributions of the carbon and oxygen stars. The ratio of the 
number of carbon stars over the number of oxygen stars 
as a function of luminosity is derived by dividing the lu- 
minosity distribution function of the carbon stars by that 
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Table 3. BVRiJK-band magnitudes, and IRAS 12 and 25 jun flux densities of the star LI-LMC1821. 

B V R i J K 5i2 (Jy) S 2S W) 

11.75 ± 0.15 10.27 ± 0.06 9.04 ± 0.04 7.52 ± 0.05 6.9 6.4 0.44 0.33 
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Fig. 9. Observed spectral energy distribution of LI-LMC1821 
(squares), and the M0 (open circles) and M5 (solid circles) 
giants from Fluks et al. (1994) scaled to a distance of 1.5 kpc, 
together with a 3000 K blackbody (dotted) 



of the oxygen stars (Fig. 10c). We have done this for the 
two extreme cases that all OC stars be carbon stars (dot- 
ted) or that all OC stars be oxygen stars (solid). The real 
luminosity relation of the fraction of carbon stars lies in 
between these two extremes. 



6. Discussion 

6.1. AGB stars 

The AGB stars that we (re-)discovered have thick CSEs 
according to their large J-K and K-[12] colours. Le 
Sidaner & Le Bertre (1996) derive an empirical relation 
between the continuum+line optical depth r of an oxygen- 
rich CSE at 10 /urn and its J-K colour: 



t - 0.8 x (J-K) - 1.2 



(9) 



Apart from the foreground star LI-LMC1821 and the 
object LI-LMC0530, the positive identifications between 
NIR and IRAS sources have J-K between 2.22 and more 
than 7.7 mag, indicating optical depths at 10 ^m in the 
range 0.6 — 5 or more. All serendipitous detections have 
J-K<2 mag, indicating small optical depths. Because we 
could not detect all stars in the J-band, the K-[12] colour 
is a better tool for studying the optical depth of the ob- 
scured AGB stars. Le Sidaner & Le Bertre (1996) derive 



Table 4. Bolometric magnitudes and object classes of the pos- 
itive identifications in the LMC: oxygen (O), carbon (C), or 
inconclusive (CO) mass-losing AGB stars, and post-AGB or 
Thermal Pulse (post/TP) stars. 



LI-LMC 


M bol 


O 


OC 


C 


post/TP 


0099 


-5.42 






X 




0109 


-6.28 


X 








0136 


-4.95 






X 




0180 


-5.01 


X 








0297 


-6.36 


X 








0326 


-5.27 








X 


0344 


-5.82 


X 








0530 


-5.55 








X 


0603 


-5.65 




X 






0777 


-4.99 








X 


0782 


-5.50 


X 








1092 


-5.57 




X 






1198 


-5.30 






X 




1316 


-4.83 








X 


1624 


-4.82 








X 


1721 


-5.27 








X 


1803 


-5.13 








X 


1813 


-5.78 




X 






1817 


-6.12 






X 




1818 


-5.89 






X 





an empirical relation between the optical depth of the CSE 
at 10 /im and the K-[12] colour: 



logr = 4.51 x log (K - [12]) - 3.37 



(10) 



Our obscured AGB stars would have optical depths be- 
tween 0.6 and 45. Post-AGB stars would not satisfy both 
the J-K and K-[12] relations for the optical depth, as 
these two relations define a path in the K-[12] versus J-K 
diagram, that is unique to the obscured AGB stars: 



(K - [12]) 4 - 51 = 1875 x (J - K) - 2813 



(11) 



It is clear from Fig. 11 that Le Sidaner & Le Bertre give a 
sequence for oxygen-rich obscured AGB stars which is rea- 
sonably consistent with our oxygen star sequence. We note 
that they corrected the IRAS 12 ^m fluxes for the spectral 
slope, assuming blackbody temperatures between 300 — 
2000 K. The colour correction factor then ranges between 
about 0.92 and 1.38 respectively. This makes their K-[12] 
colours change by —0.09 mag in the case of a 300 K black- 
body, and 0.35 mag in the case of a 2000 K blackbody. As 
the stars with thick CSEs generally have the lower black- 
body temperatures, the difference between the Le Sidaner 
& Le Bertre track and our oxygen star sequence in the 
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Fig. 10. (a) Cumulative distribution function of the stars 
over absolute bolometric magnitude: "oxygen" stars (solid), 
"oxygen or carbon" stars (dashed), "carbon" stars (dot- 
ted), all AGB stars (bold solid), and post-AGB candidates 
(dash-dotted), (b) Applying Gaussian broadening (see text), 
we derive the distribution functions N. (c) From this we derive 
the distribution function of the ratio of carbon to oxygen AGB 
stars for the two extreme cases of carbon star favoured (all OC 
stars are C stars; dotted) and oxygen star favoured (all OC 
stars are O stars; solid) 



K~[12] versus J-K diagram would only be appreciable at 
the smallest J-K, and only a few tenths of a magnitude 
at most. 



n i I r 




- carbon stars 

- oxygen stars 

I,e Sidaner & I,e Bertre 



2 4 6 8 

J-K 

Fig. 11. K-[12] versus J-K colours according to our sequences 
for carbon (dotted) and oxygen (solid) obscured AGB stars, 
and according to the parameterisation of the optical depth at 
10 /j,m by Le Sidaner & Le Bertre (1996, dashed). All are de- 
rived empirically. We adopt [12]= -2.5 log (Si 2 /28.3), with S 12 
the flux density in Jy in the IRAS 12 /im band (IRAS Explana- 
tory Supplement 1988) 



The most luminous star in our sample has a bolomet- 
ric magnitude Mboi = —6.4 mag, a factor of two fainter 
than the theoretical AGB limit as derived from the Chan- 
drasckhar limit for the core mass and the core mass- 
luminosity relation from Paczyriski (1971). The mass- 
losing AGB star sample of paper II spanned the lumi- 
nosity range between M bol <~ —4.9 and —7.2 mag. Our 
new sample consists of fainter stars on average, but with 
Mboi = — 5.0 mag our faintest star is not fainter than those 
of paper II. 

Carbon stars are more numerous at fainter luminosi- 
ties relative to oxygen stars. The two brightest stars have 
oxygen rich CSEs, but the brightest carbon star, with 
Mboi = —6.1 mag, is not much fainter than these two 
oxygen stars. On the other hand, the faintest star is a 
carbon star, but the faintest oxygen star has the same 
Mboi = —5.0 mag. Hence we conclude that we do not de- 
tect a luminosity regime in which the mass-losing AGB 
stars are exclusively either carbon stars or oxygen stars. 

This could mean that luminous AGB stars are pre- 
vented from becoming carbon stars, starting from lumi- 
nosities as low as M bo i ~ —5 mag. Hot Bottom Burning 
(HBB) may be responsible for this, as the inner boundary 
of the convective mantles of the more massive AGB stars 
becomes sufficiently hot for CN processing to occur (Sugi- 
moto 1971; Iben 1975; Scalo et al. 1975). The co-existence 
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of oxygen and carbon stars over a large range of luminosity 
may be explained as a consequence of a spread in metal- 
licity. Mass-losing oxygen stars would be metal poor and 
experience HBB at Mboi ~ — 5 mag, whereas mass- losing 
carbon stars would be more metal rich and not experience 
HBB at Mboi ~ —5 mag. A thorough discussion of the lu- 
minosity distribution function and mass-loss rates of the 
mass-losing AGB stars in the LMC is postponed to the 
next paper in this series, in combination with the sample 
of paper II. 

6.2. Future searches for mass-losing AGB stars 

We have been successful in detecting NIR counterparts 
in approximately two out of every three cases. Of these, 
approximately two out of every three cases turned out to 
be mass-losing AGB stars. Of these, approximately one 
out of every two cases was too faint to be detected in the 
J-band. Would it still be worthwhile to search for NIR 
counterparts of the remaining 37 IRAS sources? Current 
attempts to find NIR counterparts may in some cases al- 
ready be limited by the ground-based NIR searches rather 
than by the IRAS mid-infrared detections. The ISO mis- 
sion is expected to yield an extensive data base of new 
mid- infrared point sources. If ISO will detect bolomet- 
rically fainter stars exhibiting similar mass-loss rates to 
those of the obscured stars which are IRAS counterparts, 
then it will be difficult or impossible to detect their NIR 
counterparts with presently available ground-based instru- 
ments. However, if it turns out that (nearly) all of the new 
ISO detections have NIR counterparts, it will demonstrate 
that bolometrically fainter stars exhibit lower mass-loss 
rates. 

6.3. Post- AGB stars 

To investigate further the nature of the suspected post- 
AGB candidates in our sample, we have taken a closer 
look at the VHG-89 sample. They classify their post- AGB 
candidates into classes I, II, III, iva, and ivb, based on 
the shape of the infrared spectral energy distribution. The 
higher the class, the more evolved the post-AGB object: 
class I is characterised by warm dust completely obscur- 
ing the underlying cool star, while in later classes the dust 
shell becomes detached, exposing the underlying, increas- 
ingly hotter star. They also argued that classes I and II 
result from more massive stars than the later classes. For 
each class of stars in the VHG-89 sample we have calcu- 
lated the mean and standard deviation of the J-K, K- 
[12], and [12]— [25] colours. We drew boxes in the colour- 
colour diagrams, centred at the mean colours and having 
sides measuring two times the standard deviations of the 
colours. In this way we obtain a schematic picture of the 
positions of the VHG-89 classes in the colour-colour dia- 
grams (Fig. 12). Indeed, the dust becomes cooler (larger 
[12]-[25] colour) and optically thinner (smaller J-K) with 



increasing class. The early classes with obscured stars gen- 
erally have larger K-[12] and/or larger J-K than do the 
optically visible stars from class IV. 

The LMC post-AGB star candidates could be identi- 
fied with the VHG-89 classes n and in. The dust of the 
LMC stars may be warmer than that of similar stars in 
the Milky Way if the dust-to-gas ratio in CSEs of LMC 
stars is smaller, permitting stellar radiation to permeate 
farther out into the CSE, thus heating the dust to higher 
temperatures. Also the water abundance may be smaller 
in the LMC, causing the oxygen-rich CSEs to cool less 
efficiently Similarly, lower CO and/or HCN abundances 
may cause warmer carbon-rich CSEs. This would result 
in somewhat smaller [12]— [25] in the LMC than in the 
Milky Way. The colours of the LMC post-AGB star can- 
didates are indeed similar to the VHG-89 stars with rel- 
atively warm dust. The LMC post-AGB star candidates 
may have slightly smaller J-K colours because their CSEs 
are optically thinner due to the lower dust-to-gas ratio, 
but this could also be explained by selection effects: three 
of our newly identified IR stars with large K-[12] colours 
have lower limits to their J-K colours that would still per- 
mit them to be post-AGB candidates. On the other hand, 
the VHG-89 sample is constructed with a blue cut-off at 
[12]— [25] = 2 mag, selecting against the bluemost sources. 

We compared the luminosities of the LMC post-AGB 
star candidates to those of the five R Coronae Borealis 
(RCB) stars known in the LMC (Alcock et al. 1996), which 
are also believed to be post-AGB stars, and to the lu- 
minosities of PNe in the LMC. The latter we took from 
Dopita & Meatheringham (1991), Zijlstra et al. (1994), 
and Dopita et al. (1997). For some of these PNe we know 
whether they result from an oxygen- or carbon-rich AGB 
star. The luminosity distributions are presented in Fig. 13. 
We note that we do not have complete samples of any of 
the type of objects shown. The distributions of the post- 
AGB star candidates are only shown down to Mboi = —4.6 
mag, since this is approximately the detection limit of the 
combined IRAS-IRAC2b search. The distributions of the 
post-AGB stars, RCB stars, and PNe are very similar, all 
dropping steeply from Mboi ~ —5 to Mboi ~ —6 mag. In 
fact, the only object more luminous than Mboi = —5.6 
mag is the unusual PN SMP-83 (Dopita et al. 1993). 

The PNe for which the C/O ratio is known can be used 
to derive the distributions of the oxygen- and carbon-rich 
PNe in the LMC, and their ratio (Fig. 13c). From this 
we conclude that the chemical composition of the PNe is 
consistent with the chemical composition of the present 
small sample of new mass-losing AGB stars. 

We note that some of the PNe and two of the RCB 
stars in the LMC have been detected by IRAS at 12 /<m 
at a level of ~0. 1—0.2 Jy (Zijlstra et al. 1994; Moshir et al. 
1992; Alcock et al. 1996), indicating that objects are ca- 
pable of maintaining their 12 /xm flux after they have left 
the AGB. Hence the fact that the LMC post-AGB star 
candidates are detected at 12 /im by IRAS does not nee- 
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Fig. 12. K-[12] colours versus [12]-[25] (a) and J-K 
(b) colours, and J-K colours versus [12]— [25] colours (c), 
for the VHG-89 sample of post-AGB stars. We fol- 
low their classification into classes I, II, III, iva, and 
ivb. The boxes are derived from their data by calcu- 
lating the mean and standard deviation for each colour, 
per object class. We adopt [12]= -2.5 log (£i 2 /28.3) and 
[12]-[25]= -2.51og(S , i 2 /S , 25) x (6.73/28.3), with 5i 2 and S 25 
the flux density in Jy in the IRAS 12 and 25 fim bands, re- 
spectively (IRAS Explanatory Supplement 1988). The Galactic 
AGB carbon star (dotted) and oxygen star (solid and dashed) 
sequences are as in Fig. 7 and Fig. 8 



essarily imply their post-AGB age. Zijlstra et al. note that 
the IRAS detected LMC PNe have blue [12]-[25] colours, 
relative to Galactic PNe. They attribute this to selection 
effects. The LMC post-AGB star candidates also have blue 
[12]-[25] colours, relative to Galactic post-AGB star can- 
didates. NIR-optical spectroscopy for the faint, obscured 
LMC post-AGB star candidates is difficult, but may be 
possible. If they are indeed post-AGB stars, their spectra 
are expected to be of intermediate type (A, F, G). How- 
ever, this would not exclude the possibility that they are 
binary systems (e.g. Whitelock et al. 1995). 

If we assume that, down to a certain lower bolometric 
luminosity limit we are equally incomplete for the cur- 
rently available samples of mass-losing AGB stars (com- 
bining paper II with this paper), post-AGB stars, and 
PNe then we can, in principle, estimate their relative 
lifetimes. It is more difficult to do this compared to the 
AGB stars that have not been detected by IRAS, because 
many of them will still evolve significantly in bolometric 
luminosity. Comparing bolometric luminosity limited sam- 
ples of IRAS detected and non-detected AGB stars there- 
fore results in comparing different populations of stars, 



with different main-sequence masses. A synthetic evolu- 
tion approach is needed to infer the AGB lifetimes for 
stars of different main-sequence masses (Groenewegen & 
de Jong 1994). The known PNe have been selected in a 
different way than the post-AGB and mass-losing AGB 
stars, because of the different observational properties of 
these objects. Hence it is not obvious that the currently 
known sample of PNe is equally incomplete as the cur- 
rently known samples of post-AGB and mass-losing AGB 
stars. Amongst the IRAS point sources that remain to be 
searched for NIR counterparts, there may be a significant 
number of post-AGB and/or mass-losing AGB stars, in 
perhaps different relative numbers than those presently 
observed. There may also be none. Amongst the IRAS 
point sources that we could not identify with a NIR coun- 
terpart, there may be optically visible post-AGB stars, or 
not. The most we can say is that the currently available 
data suggest that the mass-losing AGB, post-AGB, and 
PN stages all have similar lifetimes, for a star with an 
AGB-tip bolometric luminosities between Mboi ~ — 5 and 
—6 mag, i.e. that has a progenitor mass in the range 2.5-4 
M Q (Vassiliadis & Wood 1993). More massive stars seem 
to have post-AGB and PN lifetimes that are considerably 
shorter than their mass-losing AGB lifetimes. 

6.4- Thermal Pulse stars 

The suspected secondary sequence in the K-[12] versus J- 
K diagram for Galactic oxygen stars does have some over- 
lap with the VHG-89 stars, but the latter are mostly found 
having larger K-[12]. This can be partly due to larger opti- 
cal depths of the CSEs of the VHG-89 stars, to the extent 
that the CSE becomes optically thick in the K-band, but 
it cannot explain the position of the VHG-89 class iv at 
small J-K. Probably the VHG-89 stars have more mas- 
sive CSEs, yielding a larger 12 /an flux, but these CSEs 
are detached, yielding small column densities towards the 
star and hence small J-K. Consequently, stars at the sec- 
ondary sequence are not expected to evolve directly into 
stars with larger K-[12]. If the stars from VHG-89 are 
post-AGB stars, then they are expected to have evolved 
from large J-K and K-[12] , first getting smaller J-K before 
also getting smaller K-[12]. Hence the secondary sequence 
may be related to the VHG-89 class 1 only. The fact that 
we noticed a secondary sequence, rather than a gradient 
from the primary sequence into the VHG-89 areas of the 
K-[12] versus J-K diagram, suggests that the secondary 
sequence and perhaps the VHG-89 class 1 are not related 
to the VHG-89 post-AGB classes 11 — iv. We can explain 
them instead as stars that only temporarily stopped losing 
mass, possibly as a result of a thermal pulse (TP). They 
are then expected to return to the AGB and resume heavy 
mass loss (Zijlstra et al. 1992). 

It is difficult to discern a TP star from a post-AGB 
star, since both experience the same phenomenon: the 
CSE becomes detached. But at least statistically there 
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Fig. 13. (a) Cumulative distribution function of the LMC ob- 
jects over absolute bolometric magnitude: post-AGB star can- 
didates (dot-dashed), R Corona Borealis stars (dashed), oxy- 
gen-rich PNe (solid), carbon-rich PNe (dotted), and all PNe 
(including the PNe for which the C/O ratio is not known; bold 
solid), (b) The distribution functions derived in the same way 
as in Fig. 10. (c) The distribution function of the ratio of car- 
bon-rich to oxygen-rich PNe 



are differences to be expected, that can be observed in the 
near- and mid-IR: the K-[12] is expected to be statistically 
larger for a post-AGB star than for a TP star. There are 
three reasons for this. First, post-AGB stars shrink while 



maintaining the same bolometric luminosity, and therefore 
they must increase in effective temperature. This results 
in increased heating of the CSE, counteracting at least 
partly the cooling of the CSE as a result of its increasing 
distance to the star. TP stars do not change as much in 
effective temperature, and hence the detaching CSE of a 
TP star cools more rapidly than that of a post-AGB star. 
Second, as the effective temperature increases but not the 
bolometric luminosity, the star will become fainter in the 
K-band, increasing its K-[12] colour. Although TP stars 
may have decreased in bolometric luminosity, this is ex- 
pected to be at most one magnitude (Vassiliadis & Wood 
1993). Third, post-AGB stars are more evolved than TP 
stars. They have experienced mass loss for a longer time 
span, yielding more massive CSEs and consequently larger 
12 /zm fluxes. 

LI-LMC0530 has too small a [12]-[25] colour for its K- 
[12] colour to be a normal post-AGB star candidate, rela- 
tive to both the VHG-89 stars and the six LMC post-AGB 
star candidates discussed above. It is also much brighter 
in the K-band than the other LMC post-AGB star can- 
didates, which suggests an optically thin CSE. Whitelock 
et al. (1995) found two stars in the South Galactic Cap 
that they explained as TP stars, because of their lack of 
variability. These two stars have J-K and [12]-[25] colours 
similar to those of LI-LMC0530, but their K-[12] colours 
are much smaller than that of LI-LMC0530, and even 
smaller than the secondary sequence in the K-[12] versus 
J-K diagram. In fact these Galactic TP star candidates 
are more similar to normal mass-losing AGB stars. This 
can be understood if they experienced one of their first 
TPs, before they had built up a massive CSE that yields 
a significant 12 fjxa excess. LI-LMC0530 seems to take 
an intermediate position between post-AGB stars and TP 
stars. 

LI-LMC0530 is identified with the LPV SHV05 10004- 
692755: an LPV with an I-band amplitude of 1.24 mag and 
a period of 169 days. Hughes & Wood (1990) give a spec- 
tral type of M6. It obeys the period-luminosity relation for 
the K-band magnitude perfectly, although bolometrically 
it is 4 — 5 times more luminous than the AGB stars that 
define the period-luminosity relation (Reid et al. 1995). 
LI-LMC0530 must have been a mass-losing AGB star in 
the past, because it has a large 12 /zm excess. It may now 
be pulsating in a (high) overtone, and maybe it will re- 
sume Mira variability and mass loss. If it is a post-AGB 
star at present, then it must be in a very early post-AGB 
stage, because its effective temperature has not increased 
much. But in this case we would expect the star to be 
obscured by an optically thick CSE. This contradicts its 
small J-K, unless we invoke a highly asphcrical shape of 
the CSE. We conclude that the nature of LI-LMC0530 is 
uncertain: it may be an early post-AGB star of class iv, or 
a star recovering from the effects of one of its last thermal 
pulses on the AGB. Alternatively, the measured flux at 
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wavelengths shorter than ~ 2/im may have been affected 
by the presence of another star in the linc-of-sight. 

7. Summary 

Our search for NIR counterparts of IRAS point sources 
in the direction of the LMC was aimed at finding ob- 
scured AGB stars in the LMC: we found 13 of them. 
The K-[12] versus J-K diagram is shown to be well suited 
for discriminating between mass-losing oxygen and carbon 
stars. Comparison between the luminosity distributions of 
the oxygen and carbon stars suggests that the fraction of 
carbon stars is smaller at bolometric luminosities around 
Mboi ~ — 6 mag than around Mboi ~ — 5 mag. Oxygen 
and carbon stars co-exist at all luminosities, but the new 
sample of obscured AGB stars does not include luminosi- 
ties as bright as Mboi ~ —7 mag. Analysis of the combined 
sample of known obscured AGB stars in the LMC will be 
presented in the next paper of this series. 

Besides the obscured AGB stars, we found a probable 
Galactic halo star or very luminous LMC red supergiant. 
We also found 7 probable IRAS counterparts which have 
large J-K colours, but not as large as the obscured AGB 
stars. Comparison with a sample of post- AGB objects in 
the Milky Way as compiled by van der Veen et al. (1989) 
suggests that these NIR detections are likely to be post- 
AGB stars, with luminosities around Mboi ~ —5 mag. We 
investigated the possibility that they are AGB stars re- 
covering from the occurrence of a thermal pulse. Although 
this cannot be ruled out completely, we found it difficult to 
reconcile their NIR and IRAS colours with those of stars 
in the Milky Way that are suspected to experience this 
short evolutionary phase. 

We also detected several objects with J-K colours 
larger than for normal stars, but not as large as the IRAS 
counterparts. They are not related to the IRAS source in 
the field. Although two of them arc identified with known 
LPVs in the LMC, the remaining few stars may be LMC 
stars that are reddened by interstellar extinction inside the 
LMC. This suggestion is supported by the serendipitous 
detection of a few galaxies with larger J-K colours than 
normal galaxies, indicating reddening by the LMC corre- 
sponding to a visual extinction of a few mag. This would 
have severe consequences for the study of stars inside the 
LMC, that are often assumed not to suffer significant in- 
terstellar extinction. 
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